Low-density expansion for the two-dimensional electron gas 
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(Dated: August 23, 2006) 

We show that in two dimensions {2D) a systematic expansion of the self-energy and the effective 
interaction of the dilute electron gas in powers of the two-body T-matrix Tq can be generated from 
the exact hierarchy of functional renormalization group equations for the one-particle irreducible 
vertices using the chemical potential as flow parameter. Due to the interference of particle-particle 
and particle- hole channels at order Tq, in 2D the ladder approximation for the self-energy is not 
reliable beyond the leading order in Tq. We also discuss two-body scattering in vacuum in arbitrary 
dimensions from the renormalization group point of view and argue that the singular interaction 
proposed by Anderson [Phys. Rev. Lett. 65, 2306 (1990)] cannot be ruled out on the basis of the 
ladder approximation. 

PACS numbers: 71.10.Ca, 71.10.Hf, 73.43.Nq 



Calculating the physical properties of strongly inter- 
acting electrons in two spatial dimensions (2D) remains 
one of the big challenges of condensed matter theory. 
Most authors agree that the normal state of interacting 
electrons in 2D is a Fermi liquid, at least for weak in- 
teractions. There are even rigorous proofs that at weak 
coupling certain two-dimensional models for interacting 
electrons (including the repulsive Hubbard model away 
from half filling) are Fermi liquids above an energy scale 
which is non-perturbative in the coupling constant 
However, due to a lack of controlled methods, one cannot 
exclude the possibility that for sufficiently strong interac- 
tions Fermi liquid theory breaks down in 2D. Anderson 
has argued repeatedly [3, S| that due to non-perturbative 
effects which are neglected by the usual field-theoretical 
machinery of many-body theory, the normal state of the 
two-dimensional Hubbard model is not a Fermi liquid 
for arbitrary strength of the interaction. This scenario 
has been criticized early on by Engelbrecht and Rande- 
ria 0], who pointed out that a calculation of the self- 
energy E(fc,a-') within the ladder approximation (LAP) 
predicts Fermi liquid behavior. However, the LAP has 
some rather peculiar and most likely unphysical features; 
in particular, the limiting behavior of Y.{k,uj) for small 
frequencies uj and for wave- vectors k in the vicinity of the 
Fermi surface kp is very sensitive to the order in which 
the limits lu ^ and k kp are taken fs, 6]. 

Unfortunately, going beyond the LAP is very difficult, 
because the particle-hole scattering channels have to be 
taken into account on the same footing with the particle- 
particle scattering channel retained within the LAP, i.e., 
one has to solve coupled Bethe-Salpeter equations in sev- 
eral channels. In this work we shall take a fresh look at 
this problem using a formally exact hierarchy for renor- 
malization group (RG) flow equations for the one-particle 
irreducible vertices |7|. 



Functional RG in the spin-singlet channel. Start- 
ing point of our analysis are the formally exact RG 
flow equations for the irreducible self-energy 'Sa{K) 
and the antisymmetrized effective interaction vertex 
r'-^\K[alK'2a'^;K2a2,Kiai). Here K = (fc,iw) de- 
notes wave- vector and Matsubara frequency, and Ci are 
the spin projections. A is a cutoff scale separating the 
low-energy from the high-energy degrees of freedom, such 
that for A = Ao the low-energy fluctuations are sup- 
pressed. For the Hubbard model on a D-dimensional 
hyper-cubic lattice with lattice spacing a and on-site re- 
pulsion U the bare interaction is F^^-* = {5ai.(j{5a2.<j'^ ^ 
'5o-i.(T^'5o-2.crj)C/o, where Uq = a^U. As we iterate the RG, 
the Ifowing effective interaction develops a spin-triplet 
component, so that in general 



+ i6,,,.'^6^,,.'^+S^^yJ,^^^,jT'^iK[,K:,;K2,K,) .(1) 

The singlet part r^(_ft'(, ifj! ^2, Ki) is symmetric under 
the exchange K[ ^ K'2 or Ki ^ K2, while the triplet 
part r'j^{K[, K2, Ki) is antisymmetric and therefore 
vanishes for K[ — K!^ or Ki — K2. The RG flow equa- 
tion for rj^'' contains also mixed terms F^F^ and a term 
involving the six-point vertex F^-*. However, F^ and F^'' 
are both irrelevant in the RG sense 0], so that we shall 
neglect them. Within this approximation, the RG flow 
of the self-energy is determined by 



dtj:„{K)^ GA{K')mK,K';K',K) , (2) 
Jk' 

and the effective interaction satisfies at zero temperature 
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dArXiK[,K:,;K2,Ki) = - [ [Ga{K)Ga{Ki+K2-K) + Ga{K)Ga{Ki+K2-K)] 

J K 

xTl{K[,K'2;Ki + K2-K, K)TX{K, Ki + K2 - K; K2, Ki) 
- \ J^{[Ga{K)Ga{K + Ki-K[) + Ga{K)Ga{K + K, - K[)] 

xTX{K[,K + K,-K[;K,Ki)TI{K:2,K;K + Ki-K[,K2) + {K[^K'2)} , (3) 



K J (27r)" 27r ' 



and we have introduced the 



where 

cutoff-dependent propagator 

GAiK) = &Aimiu - a - ^AiK)Y 



(4) 



and the single scale propagator 

Ga{K) = [1 + Gl{K)J:A{K)]-^dAGl{K) . (5) 

Here 0a (fe) is some cutoff function which suppresses the 
low-energy modes, normalized such that Oo(fc) = 1 and 
©ooCs) = 0. The cutoff-dependent non-interacting prop- 
agator is G^K) = QA{k)[iuj - ffc]"\ with = efc - /i, 
where Ck is the bare energy dispersion and /i is the chem- 
ical potential. The first term on the right-hand side 
of Eq. ^ is the contribution from the particle-particle 
channel, while the other two terms correspond to the 
particle- hole channels (also called zero-sound channels). 

Ladder approximation. In the limit of vanishing den- 
sity only the particle-particle channel contributes to the 
effective interaction. At low densities, it is then rea- 
sonable to calculate the self-energy in LAP, where the 
contribution of the particle-hole channels to the effective 
interaction in Eq. ^ is simply ignored. It is then advan- 
tageous to introduce total and relative energy-momenta 
as independent variables, P — Ki + K2 = K[ + K2, 
Q = {Ki- K2)/2, Q' = {K[ - K'^)/2, and define 

rA(|+Q',f +Q) = rA(Q',Q;P) . (6) 

With this approximation Eq. fS]) reduces to 

dATA{Q',Q;P) - -2 

xTa{Q\K-P)Ta{K,Q-P) . (7) 



Ga{^+K)Ga{^-K) 



Assuming that the bare interaction is independent of 
the frequency-part of Q and Q', this remains also true 
for the renormalized interaction, so that we may write 
Ta{Q\Q;P) = rA{q',q;P)- If we replace the flowing 
Green functions in Eq. ([7]) by the non-interacting ones, 
the frequency sum is easily carried out. 



dATA{q',q;P) ^ j ^ 
I- f{U+k)- f{U-k) 



^0aA[eA(f + fc)eA(f-fc)] 



+fc - ^f-* 



TA{q\k-P)VA{k,q-P) , (8) 



where P — (p, iw), and /(^) = 9 
tion at zero temperature. Eq. ([5 
Bethe-Salpeter equation for the 
the particle-particle channel [8]. 
if we assume that rA{q\q;P) — 
of q and q' , in which case Eq. ([B 
differential equation, dA^A{P) — Xa{P)^a{P) , where 

-iDi 



(— ^) is the Fermi func- 
I is equivalent with the 
effective interaction in 
We recover the LAP 
Ta{P) is independent 
becomes an ordinary 



xa(^) = / ^aA[eA(| + fc)eA(|-fc)] 



iu — £v 



+fc 



(9) 



The solution of the above differential equation with ini- 
tial condition Taq (P) — Ua yields the usual LAP for the 
effective interaction, 

r(p) = TaMp) = Uo[i + Uox{P)r' , (10) 

where the pair susceptibility x{P) = Jq" d,AxA{P) is 
given by 



xiP) 



d^k 



i-eA„(f + fc)eA„(|-fc) 



(11) 



+k 



Zero density limit. In this limit /i ^ and we may 
approximate Ck = fc^/(2m), where m is some effective 
band mass. The Fermi functions /(Cfifc) then vanish. 
Using for simplicity a sharp cutoff in momentum space, 
eA(fe) = e(|fc| - A), we find from Eq. © that the two- 
body T-matrix TA{q' , q; iCo) = rA(q', q; p = 0, iui) in vac- 
uum satisfies the fiow equation 



dATA{q',q,iuj) 



d^k <5(|fc|-A) 
(27r)^ zw- A2/TO 
xTA{q\k,iQ)TA{k,q,iLo) . (12) 



Uq'-q, where 



The initial condition is TA„iq' , q,iLo) 
Uk is the Fourier transform of the bare interaction. 
Eq. (fT2|) is the RG version of the Lippmann-Schwinger 
equation for the T-matrix of elementary scattering the- 
ory [9]. The low-energy behavior of the T-matrix de- 
fines the s-wave scattering length via r(0, 0, iO) = 
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(Color online) RG flow of the rescaled zero energy 
= mKDA°''^TA{0, 0, iO) of the two-body T-matrix. 



FIG. 1: 

part ui — 

7_DCif~^/m, where 7^) is a numerical constant (73 = 
47r). Introducing the dimensionless coupling constant 
ui = mKDA^~^TA{0,0,iO) = Ko-fDiasA^)"-^ , where 
I = ln(Ao/A) and Kd = /T{D /2) is the sur- 

face area of the D-dimensional unit sphere divided by 
{2tt)^ , we obtain from Eq. the RG flow equation 



dm = eui 



2- D 



(13) 



The corresponding flow of ui is shown in Fig. [T] For 
D > 2 the RG flow has a stable Gaussian fixed point 
G at vanishing interaction, and an unstable fixed point 
UF at finite negative interaction = e < 0. At this 
point the zero energy T-matrix diverges for A — > 0, 
Ta(0,0,iO) oc u*A^~^ —00, so that UF corresponds 
to the unitary Fermi gas in the limit of vanishing den- 
sity For D — 2 both fixed points UF and G merge 
into a single point UFG at zero interaction. Finally, for 
D < 2 the unitary fixed point UF has zero interaction, 
because in this case an infinitesimal attractive interaction 
leads to a zero energy bound state fvi^. However, there is 
a new non-Gaussian fixed point NG for finite repulsive in- 
teraction = e > 0, which resembles the Wilson- Fisher 
fixed point below four dimensions We now show 

that this fixed point is characterized by a finite phase 
shift 6o{0) — — 7re/2 for s-wave scattering at zero energy. 
Setting again Uk = Uq, the solution of Eq. is 

To(icS) - TA^oiq',q,iOj) - Uo[l + UoXo{ii^T^ , (14) 

where from Eq. (jlip we obtain 

d^k e(Ao - |fc|) 



(2n)D iuj-k^lm 



(15) 



For D < 2 the integral is convergent for Aq ^ cxd, and 
we obtain for real frequencies, 

Xo(w + «0) = 7rt/£)(|a;|, r7ir){9(— cj)[sin(7re/2)]^"^ 

-f-e(w) [cot(7re/2) -f ^]} , (16) 

where vj^iuj^nir) — K]jmj.{2mj.uj)^'^/'^ is the density of 
states of a free particle with reduced mass — m/2. 
In 2D, where the density of states V2{m) ~ m/(27r) is 
constant, we have 



where Eq = Aq/tti. The phase shift for s-wave scattering 
in the spin-singlet channel, defined via [i, 13| To(a;-|-iO) — 
-|To(u; -I- iO)|e'*n("), is for w > given by 



(5o(w) 



ar 



ctan{7rg(tj)/[l + 7rg(tj)cot(7re/2)]} , (18) 



where g{Lo) = V£i{u)t 711^)110 is the relevant dimension- 
less coupling. Because .g(0) = cxd for Z? < 2 due to 
the divergence of the density of states, we conclude that 
(5o(0) = — 7re/2 for any finite Uq > 0. 

Low-density expansion in 2D. If we approximate the 
effective interaction in Eq. ^ by r(P) given in Eq. (fTU]) 
and ignore self-energy corrections to the propagators, we 
obtain the usual LAP for the self-energy [J, l5|. However, 
the LAP is not a systematic low-density expansion. In 2D 
the relevant dimensionless low-density parameter is a = 
[ln(Ao/fcF)]~"'^. One should therefore expand the LAP in 
powers of a. The self-energy within LAP up to order 
has been evaluated by Bloom [iJl using the pseudo- 
potential method pioneered by Galitskii [1[ . In this 
method one first derives a non-linear integral equation 
relating the effective interaction in LAP to the on-shell 
two-body T-matrix in vacuum. Iteration of this integral 
equation generates the desired expansion in powers of a. 

The RG equations ([2]) and ([3]) offer an alternative way 
to generate the low-density expansion in 2D. The ad- 
vantage of the RG approach is that the contribution 
from the particle-hole channels neglected in LAP can 
be systematically included. The crucial point is that 
Eqs. ([2]) and ^ remain formally valid if we take the 
derivative with respect to any parameter appearing only 
in the free propagator Q. In particular, we may choose 
A = /i, in which case Eq. ^ and ^ describe the 
evolution of the self-energy and the effective interac- 
tion as we change the chemical potential. The propa- 
gator in the flow equations should then be replaced by 
G^{K) = [iuj — Ck + ^ — E^(iir)]~^, and the single scale 
propagator by G^{K) — ~G'j^{K). The flow equation for 
the self-energy now reads 



d^E,iK)= I G,{P-K)r,i^~Kr--K;P), (19) 



P 
"2 



and the flow equation for the effective interaction can be 
obtained by making analogous substitutions in Eq. 
Assuming that an expansion in powers of fi is possible, we 
may generate this expansion by successive differentiation 
of Eq. HH) and the analogue of Eq. For Ok = Uo 
the initial conditions at fi = are S^^o(-f^) = aud 
r,,.o(f - K, ^ -K;P)= Tq{ilo - ^ T{P). The 
leading term in the expansion of the self-energy is 



S^(if ) = Go(P - K)T{P) + 0{^i') , (20) 



Xo(w -I- iO) = 7ri'2(»7ir) [In (£'o/|'^|) + «0(w)] , (17) and from Eqs. ([3]) and (O we obtain for the effective 
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interaction 

r^{Q\Q;P) = T{P) - tiT\P)[d^x{P)V=o 
J^{[Go{K)GoiK + Q + Q') 

+Go(K)Go{K + Q + Q')]T{K + Q + |) 

xT{K + Q' + ^) + (g' ^ -Q')} + 0{^i^), (21) 

where the pair susceptibihty x{P) finite /i is given 
in Eq. (jlip . The second term on the right-hand side 
of Eq. pT|) is due to the /i-dependence of the effective 
interaction in LAP, while the last term is the contribution 
from the particlc-hole channels neglected in LAP. 

Due to the non-analytic /i-dependence of the density 
of states for D ^ 2 the above expansion in powers of /i 
is only possible in D = 2. To see this more clearly, let us 
explicitly evaluate the first term on the right-hand side 
of Eq. (pOI . The result can be written as 



Am ' 



(22) 



where Po(m) = iJLV2{m). In fact, Eq. (122]) is the leading 
term in the expansion in powers of Tq for arbitrary D 
if we identify po{iJ.) = {2/ D)p.V£)[p,m) oc /x^/^ with the 
density (per spin projection) in the absence of interac- 
tions. The non-analyticity in /i for D ^ 2 is obvious. 

If we approximate Tq by the bare interaction J7o, then 
Eq. (|22p reduces to the Fock correction to the self- 
energy. In contrast to the usual LAP 0, Q , the particle- 
particle susceptibility is evaluated at vanishing density 
in Eq. (P^ . It is instructive to explicitly calculate the 
corresponding single-particle Green function G{k, ioj) for 
D <2. The relation between fx and Fermi momentum kp 

4m' - Expanding G~^{k,iuj) 
kp we obtain the Fermi liquid form 



is then po(/i)7'o(-4^ 
for small m and I k 



G{k, iuj) « Z[iuj ~ {kF/m^){\k\ — kp)y 



with quasi-particle residue Z = 1 - 
fective mass renormalization m/m, 



(23) 



and ef- 



-9\ 



Here g = i^£i(/x, mr)?7o, and in 2D we should replace 
e — > a = [ln(Ao/fc_F)]~^. For 5 — > cx) we obtain Z ^ 1 — 
and m/m^ = 1 — £^J2. This agrees qualitatively with the 
results by Bloom [1J|, who obtained different numerical 
coefficients on front of the correction terms, which is due 
to the fact that his approximation amounts to retaining 
also the second term on the right-hand side of Eq. pTjl . 
However, to order Tg the particle-hole contributions to 
Eq. ([2T|) compete with the particle-particle channel re- 
tained in Ref. IJ], so that the approximation employed 
by Bloom is inconsistent. In this point we agree with 
Anderson Q, but for different reasons: there is no math- 
ematical mistake in Bloom's analysis, but in 2D the LAP 
is inconsistent beyond the leading order in Tq. 



In summary, we have reconsidered the low-density ex- 
pansion for the electron gas in dimensions D < 2 us- 
ing functional RG methods. At vanishing density, we 
have rewritten the Lippmann-Schwinger equation for the 
two-body T-matrix in vacuum as a RG flow equation 
and have pointed out that for D < 2 this equation has 
a non-Gaussian fixed point characterized by a finite s- 
wave phase shift. At low densities in 2D a systematic 
expansion of the self-energy and of the effective interac- 
tion in powers of the two-body T-matrix can be gener- 
ated with the help of the functional RG flow equations 
for the irreducible vertices using the chemical potential as 
flow parameter. We have argued that in 2D the LAP is 
inconsistent beyond the leading order in this expansion 
because already at order Tg the contribution from the 
particle-hole channels competes with the particle-particle 
channel and gives rise to a complicated momentum- and 
frequency dependence of the effective interaction which 
still has to be explored. Because the LAP is not reliable 
in 2D, one cannot use the LAP to rule out that the true 
effective interaction in 2D exhibits a singular dependence 
on the relative momenta, as proposed by Anderson 0, 0] . 

We thank M. Salmhofer for his comments on rigorous 
results for two-dimensional Fermi systems. 
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